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ABSTRACT 

The recently constructed Hubble diagram using a combined sample of SNLS 
and SDSS-II Type la SNe, and an application of the Alcock-Paezynski (AP) test 
using model-independent Baryon Acoustic Oscillation data, have suggested that 
the principal constraint underlying the cosmic expansion is the total equation-of- 
state of the cosmic fluid, rather than that of its dark energy. These studies have 
focused on the critical redshift range (0 < z < 2) within which the transition 
from decelerated to accelerated expansion is thought to have occurred, and they 
suggest that the cosmic fluid has zero active mass, consistent with a constant 
expansion rate. The evident impact of this conclusion on cosmological theory 
calls for an independent confirmation. In this paper, we carry out this crucial one- 
on-one comparison between the R\ x = ct Universe (an FRW cosmology with zero 
active mass) and wCDM/ACDM, using the latest high-z measurements of H(z). 
Whereas the Type la SNe yield the integrated luminosity distance, while the 
AP diagnostic tests the geometry of the Universe, the Hubble parameter directly 
samples the expansion rate itself. We find that the mo del-independent cosmic 
chronometer data prefer A h = ct over wCDM/ACDM with a BIC likelihood 
of ~ 95% versus only ~ 5%, in strong support of the earlier SNela and AP 
results. This contrasts with a recent analysis of H(z) data based solely on BAO 
measurements which, however, strongly depend on the assumed cosmology. We 
discuss why the latter approach is inappropriate for model comparisons, and 
emphasize again the need for truly model-independent observations to be used 
in cosmological tests. 

Subject headings: cosmology: cosmological parameters - cosmology: distance scale - 
cosmology: observations - cosmology: theory - galaxies 
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1. Introduction 

Modern cosmology is based on the Friedmann-Robertson-Walker (FRW) metric to 
describe the spacetime expansion of the cosmic fluid, often assumed to be shear free (i.e., 
‘perfect’), and comprised of at least three components. Two of these, matter (p m ) and 
radiation (p r ), are readily observed, while the third, a poorly understood ‘dark energy’ (pde), 
is inferred from the analysis of the distance-redshift relationship in Type la SNe (Perlmutter 
et ah 1998; Riess et al. 1998; Schmidt et al. 1998). The standard model of cosmology, 
which we refer to as ACDM when pd e is a cosmological constant with equation-of-state 
Wde = Pde/Pde = — 1, and tcCDM when it is not, adopts these elements, but the poorly 
known w^ e burdens the theory with several unknown parameters that need to be optimized 
while fitting the data. 

In spite of this limitation, tcCDM/ACDM has done remarkably well in accounting 
for the observations, though the ever improving precision of the various measurements 
is beginning to uncover tension—in some cases, actual inconsistencies—with the model 
predictions. In fact, the most recent tests of the standard model have revealed that the 
principal constraint underlying the cosmic expansion appears to be the total equation-of- 
state for the cosmic fluid rather than that of its dark energy. The observational evidence 
now favors a cosmic fluid with zero active mass, i.e., p + 3p = 0, where p = p m + p r + pd e 
and p = p m + p r + Pde are, respectively, its total energy density and pressure (see also Mclia 
2015b). This result, however, would exclude a cosmological constant as a representation of 
dark energy. 

It is therefore critically important to seriously examine the zero active mass condition 
with the highest precision permitted by current observations. In the literature, the FRW 
cosmology with zero active mass is known as the R h = ct Universe (Melia 2007; Melia and 
Shevchuk 2012). This model has been compared with wCDM/ACDM (an FRW cosmology 
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lacking this condition) in many one-on-one tests, at both low and high redshifts (see, e.g., 
Melia and Maier 2013; Wei et al. 2013, 2014a, 2014b, 2015a; Melia et al. 2015a). Thus 
far, R h = ct has been favored by the data, with model selection tools typically yielding 
likelihoods of ~ 90% versus only ~ 10% for tcCDM/ACDM. ffowever, these results have 
not yet been universally accepted, and several counterclaims have emerged in recent 
years. These may be loosely grouped together into four primary categories: (1) that R h 
(representing the gravitational radius and, therefore, also the ffubblc radius) does not really 
have any physical meaning or bearing on the observations; this claim has been made chiefly 
by van Oirschot et al. (2010), Bilicki & Seikcl (2012), and Lewis & van Oirschot (2012); 
(2) that the zero active mass condition p + 3p = 0 cannot be made consistent with the 
actual constitutents in the cosmic fluid (Lewis 2013); (3) that the measurements of H(z ) 
as a function of redshift (the primary topic of the present paper) favor the concordance 
model over = ct (Bilicki & Seikel 2012; Shafer 2015); and (4) that the analysis of Type 
la SNe also favor the concordance model over R\ l = ct (Bilicki & Seikel 2012; Shafer 2015). 
These papers, and those published in response to them (see, e.g., Melia & Maier 2013; 
Bikwa et al. 2012; Melia & Shevchuk 2012; Melia 2012, 2014, 2015a), have generated a 
very important discussion that we aim to continue here. Specifically, in § 3 below, we will 
describe at length why the choice of truly model-independent data, and their analysis using 
sound statistical practices, is of utmost importance to any serious attempt at comparing 
different cosmological models in an unbiased fashion. 

Of the previous tests favoring R h = ct over tcCDM/ACDM, two are particularly 
noteworthy. In the first of these, a combined sample of 613 supernova events from 
SNLS (Guy et al. 2010) and the newly released SDSS-II (Sako et al. 2015) was used to show 
that in an unbiased pairwise comparison, the Bayes Information Criterion (BIG) favors 
the Rh = ct LIniverse with a likelihood of « 88% versus ~ 12% for the standard model 
(Wei et al. 2015b; Melia et al. 2015b). This combined sample spans the critical redshift 
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range (0 < z < 1.2) in which the presumed transition from deceleration to acceleration is 
thought to have occurred. This outcome reaffirms the influence of dark energy, but calls 
into question the widely held belief that the Universe is currently accelerating. In other 
words, when a neutral calibration of the SN data is used to compare models, the evidence 
does not appear to support the interpretation of dark energy as a cosmological constant. 

In a second test, an application of the Alcock-Paczynski (AP) diagnostic using the 
model-independent anisotropic distribution of Baryon Acoustic Oscillation (BAO) peaks 
at average redshifts (z) = 0.57 and (z) = 2.34 has shown that the BAO data disfavor 
taCDM/ACDM at better than 2.7 a (Melia and Lopez-Corredoira 2015). In the context of 
expanding FRW cosmologies, these measurements instead strongly prefer the zero active 
mass equation-of-state, confirming the conclusion drawn earlier from the Type la SN 
Hubble diagram. 

Each of these results is quite significant in its own right; together, they make quite 
a compelling case. But given their evident impact, they need to be confirmed using 
another independent diagnostic in the critical redshift range sampled by the SNela and 
BAO observations. In this paper, we focus on several new measurements of the Hubble 
parameter H(z) within the redshift range 0 < z < 2.5, using luminous red galaxies as 
cosmic chronometers (Moresco 2015). These measurements were made without priors, 
and are therefore model independent—ideally suited to a one-on-one comparison between 
taCDM/ACDM and R h = ct. 

With the latest extension of the H(z) coverage to z ~ 2, the cosmic chronometer data 
are now fully compatible with the SNela and BAO measurements. But whereas the SN 
comparison is based on the (integrated) luminosity distance, and the AP diagnostic tests the 
predicted geometry, the H(z) measurements directly probe the expansion rate itself. Thus, 
although they sample the same redshift range, these three diagnostics test the cosmologies 
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in distinctly different ways, and each therefore contributes to the self-consistency of the 
tests in its own unique way. In this paper, we carry out the critical one-on-one comparison 
between = ct and tcCDM/ACDM using the latest cosmic chronometer data. In so doing, 
we demonstrate that this independent test also strongly favors R h = ct over tcCDM/ACDM, 
thereby strengthening the argument in favor of the zero active mass condition. 


2. The Cosmic Chronometers 


Since the Hubble rate depends on the expansion parameter a(t) according to 

H(z) = - , (1) 

a 

and (1 + z) = a 0 /a, where a 0 is the expansion parameter today, H(z) may be measured 
directly from the time-redshift derivative dt/dz using 


H(z) 


1 dz 
1 + z dt 


( 2 ) 


Measuring H(z) using the differential age of the Universe therefore circumvents the 
limitations associated with the use if integrated histories. 


Galaxies evolving passively on a time scale much longer than their age difference are 
among the best cosmic chronometers one may use for this purpose (Jimenez and Loeb 
2002). Less than 1 per cent of the stellar mass in the most massive galaxies formed at z < 1 
(Dunlop et al. 1996; Spinrad et al. 1997; Cowie et ah 1999; Heavens et ah 2004; Thomas 
et ah 2005; Panter et ah 2007). In fact, star formation ceased by redshift z ~ 3 (Thomas 
et ah 2005) in galaxy clusters, while quite generally, all systems with stellar mass over 
~ 5 x 10 11 M q ended their star formation activity by z ~ 2 (Treu et al. 2005). 

It is therefore reasonable to assume that galaxies in the highest density regions of 
clusters have been evolving passively since z ~ 2, tracing the so-called ‘red envelope,’ which 
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Fig. 1.— Twenty five (model-independent) measurements of H(z), with error bars, and a 
comparison with two theoretical models: (dashed) the optimized tcCDM cosmology, with 
best-fit parameters H 0 = 68.0 ± 7.9 km s -1 Mpc -1 , h2 m = 0.31 ± 0.16, and Wd e = — 1 O.OI^q;^; 
and (solid) the = ct Universe, with its sole parameter Hq = 63.3 ± 7.7 km s -1 Mpc -1 . 
The reduced x( lof (24 degrees of freedom) for R h = ct is 0.56. The corresponding value for 
wCDM (22 degrees of freedom) is x^ of = 0.58. The BIG favors Rh = ct over wCDM with a 
likelihood of ~ 94.3% versus ~ 5.7%. 

























samples the oldest stars in the Universe at every redshift. For this reason, these structures 
have been used as reliable cosmic chronometers in several extended studies (Stern et al. 
2010a,b; Moresco et ah 2012a,b), culminating with the most recent measurements at z ~ 2 
(Moresco 2015). 

The data shown in figure 1 were assembled from the compilations of Simon et al. 
(2005), Stern et al. (2010), Moresco et al. (2012a) and Moresco (2015). We emphasize the 
fact that all of these measurements are mo del-independent. Other kinds of measurement of 
H(z), based on the identification of BAO and the Alcock-Paczynski distortion from galaxy 
clustering, depend on how ‘standard rulers’ evolve with redshift, rather than how cosmic 
time changes with z. Unfortunately, the values of H(z) measured in these different ways are 
sometimes combined to produce an overall H(z) versus z diagram, but such compilations 
cannot be used to test different models because the second approach necessarily adopts 
a particular cosmology and is therefore model-dependent (Blake et al. 2012). Only the 
cosmic chronometer measurements are truly mo del-independent and therefore suitable for 
comparing tcCDM/ACDM to other models, such as Ry, = ct. 


3. Discussion 

In the 7?^ = ct Universe, a(t) oc t, so the Hubble parameter scales very simply with 
redshift: 

H(z) = H 0 (l + z). (3) 

Since H 0 is the sole free parameter in this model, the corresponding theoretical curve in 
figure 1 can be adjusted only by sliding it vertically. There is no freedom to alter its 
gradient with redshift. Because the redshift coverage with these measurements now extends 
to z ~ 2, this particular test is therefore very constraining on models, such as Ry, = ct, that 
do not predict a transition from deceleration to acceleration somewhere near z ~ 1. 
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Table 1: Optimized fits of H(z) 


Model 

Optimized parameters (conventional units) 

^dof 

Likelihoods 





"5 

II 

H 0 

= 63.3 ± 7.7 


0.566 

AIC: 82.9% 

KIC: 

93.0% 

BIC: 

94.3% 

uiCDM 

H 0 

= 68.0 ± 7.9, 

= 0.31 ± 0.16, w de = -0.91-°'^ 

0.579 

AIC: 17.1% 

KIC: 

7.0% 

BIC: 

5.7% 

ACDM 

Ho 

= 73.3 ± 3.5, 

= 0.28 ± 0.04, n de = 0.60 ± 0.13 

0.580 

AIC: 17.1% 

KIC: 

7.0% 

BIC: 

5.7% 


wCDM/ACDM is characterized by a larger number of free parameters because the 
properties of dark energy are so poorly known. At the very least, one must include the 
Hubble constant H 0 , the fractional energy density Q m in the form of matter (= p m [h)]/Pc| 
where p c = 3 c 2 Hq/8tiG is today’s critical density), the dark energy equation-of-state 
parameter, w de = p de /pde, and the spatial curvature constant k. But to keep the comparison 
as simple and favorable to icCDM/ACDM as possible, we will here adopt the prior value 
k — 0 in the case of icCDM, and w^e = —1 for ACDM, thus reducing the total number of 
adjustable parameters to an essential three in each case. The Hubble parameter in this 
cosmology is therefore given by the expression 

H(z) = H 0 [fi m (l + £) 3 + fi r (l + zf + D de (l + z) 3 ( 1+tt,de) ] 1/2 . (4) 

Here, h2 r for radiation and fl de for dark energy are defined analogously to Q m . The scaling 
of each density with redshift assumes that all three components evolve independently of 
each other, which may be reasonable in the local Universe (i.e., at least out to z ~ 2). 
Note that although we formally include Q r in this expression, practically speaking its value 
(~ 6 x 1CU 5 ) remains too insignificant compared to the other densities for it to contribute 
to the evolution in H{z) in this redshift range. 

For each model, we optimize the fit by finding the set of parameters that minimize the 
X 2 , Ho in the case of R h = ct; H 0 , H m and w de in the case of wCDM; and H 0 , fl m and H de 
for ACDM, as is evident from Equations (3) and (4). The results are summarized in Table 
1, and the corresponding best fits for R ^ = ct and icCDM are shown in figure 1. These 
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optimized parameter values are quoted with one-sigma standard errors, calculated from the 
corresponding ^-distribution for each model. With 25 — 1 = 24 degrees of freedom, the 
reduced Xdof f° r the -^h = ct model is 0.566. By comparison, the optimal tcCDM fit has 
25 — 3 = 22 degrees of freedom, and a corresponding x\ 0 { = 0.579. The best-fit ACDM 
model is similar to this, with \ 2 doi = 0.58. These fits suggest that Ah = ct is at least as good 
as wCDM/ACDM in accounting for the model-independent H(z) measurements, especially 
since it has only one free parameter. In particular, note how well it accounts for the most 
recent z ~ 2 measurement; though, to be fair, the error bar at this redshift is still quite 
large, so that tcCDM/ACDM cannot be excluded simply on the basis of this value alone. 

On statistical grounds, however, the A h = ct expansion rate versus redshift is far 
more likely to be correct than that in tcCDM/ACDM. To compare the evidence for and 
against competing models, the use of various information criteria is now fairly common in 
cosmology (see, e.g., Takeuchi 2000; Liddle 2004, 2007; Tan and Biswas 2012). For example, 
the Akaike Information Criterion (AIC), 

AIC = x 2 + 2/ , (5) 

provides an enhanced ‘goodness of fit,’ extending the usual y 2 statistic by also taking 
into account the number of free parameters / in each model. The AIC ranks two or 
more competing models, and yields a numerical measure of confidence that each model 
is preferred (Akaike 1973; Burnham and Anderson 2002, 2004). Clearly, the AIC prefers 
models with fewer parameters, as long as the others do not provide a substantially better 
fit to the data. Informally, in a one-on-one comparison between two models, model a has a 
likelihood 

T _ exp(—AICq/2) _ 

“ exp(-AIC 1 /2)+exp(-AIC 2 /2) U 

of being the best choice. For the analysis reported here, these likelihoods are 82.9% for 

i?h = ct and 17.1% for either wCDM or ACDM (see Table 1). 
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An alternative to the AIC (Cavanaugh 1999) is known as the Kullback Information 
Criterion (KIC), KIC = y 2 + 3/, which disfavors overfitting more than does the AIC. Not 
surprisingly, the corresponding values of KIC = 93.0% for R h = ct and 7.0% for either 
tcCDM or ACDM, therefore favor the former by a greater amount. 

A better-known alternative, known as the Bayes Information Criterion (BIC), is 
actually not based on information theory, but rather on an asymptotic approximation to the 
results of a conventional Bayesian inference procedure for comparing two models (Schwarz 
1978): 

BIC = y 2 + /(In n) . (7) 

The BIC suppresses overfitting very strongly when the number of data points n is large. In 
the sample used here, n = 25, and we find that BIC = 94.3% for R h = ct, versus only 5.7% 
for either wCDM or ACDM. 

According to all three statistics, the expansion rate versus redshift predicted by i?h = ct 
is more likely than that of icCDM/ACDM to be closer to the ‘correct’ cosmology. This 
consistency in the results is quite important because, whereas the BIC is based on Bayesian 
statistics, the AIC and KIC are not. When n is large, any chosen Bayesian ‘priors’ over the 
parameters of the individual models drop out (Kuha 2004). However, the prior distribution 
over the choice of models does not. In other words, though the model distribution is often 
reasonably assumed to be ‘flat’, i.e., to give each model an equal a priori likelihood, a case 
can sometimes be made that one model is to be preferred a priori over the other. 

The idea is that in the light of the data, e.g., the BIC computed from the data, the 
relative likelihood of the two models being compared is multiplied by a Bayes factor. When 
one assumes that the models are equally likely a priori, the multiplication by the Bayes 
factor yields (if the BIC is being used) their likelihoods being given by the BIC-derived 
Bayes weights, according to the usual formula. But the models don’t necessarily need to be 
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taken, a priori (i.e., pre-data), to be equally likely. 

So for the BIC one needs to assess beforehand whether either taCDM/ACDM or 
i?h = ct is preferred. Of course, this is precisely the problem with attempting to prioritize 
models before testing them against the data. Different workers can have different subjective 
opinions. As of today, these two models have been tested against each other using diverse 
types of data, all of which have thus far favored = ct. Should one assign a prior 
distribution for the models, the evidence now suggests that Rh = ct is preferred. But to 
avoid such biases based on subjective points of view, the only sensible approach with the 
work reported in this paper is to simply let the data speak for themselves, and assume that 
the prior distribution over the choice of models is flat. 

The fact that we here use three different information criteria to evaluate the models, of 
which the AIC and KIC have their own foundations and do not belong to Bayesian statistics, 
strengthens our overall conclusions. The inclusion of new liigh-z measurements of H(z) 
using cosmic chronometers has reinforced the results of cosmological model comparisons, all 
of which have thus far favored R h = ct over tcCDM/ACDM. 

Nonetheless, in a similar study of H(z) measurements based on BAO observations, 
rather than cosmic chronometers, Shafer (2015) reaches a different conclusion. Unfortunately, 
Shafer’s analysis is a good example of the unwitting use of mo del-dependent measurements 
to test competing cosmologies. We reiterate the obvious, though often forgotten, necessity 
of using only truly model-independent data to conduct statistically fair comparisons 
between different models. Without adequate explanation, Shafer opts to completely ignore 
measurements of H(z) using cosmic chronometers, relying instead on BAO measurements, 
in spite of the fact that the deficiencies in the latter have already been discussed at length, 
e.g., in Melia & Maier (2013). He appears to be unaware of the significant limitations of all 
but the most recent two BAO measurements at z = 0.57 (Anderson et al. 2014) and 2.34 
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(Delubac et al. 2015) for this kind of work. All previous applications of the galaxy two-point 
correlation function to measure a redshift-dependent scale were limited by the difficulty 
in disentangling the acoustic length in redshift space from redshift distortions due to 
internal gravitational effects (Lopez-Corredoira 2014). With this process, one had to either 
pre-assume a particular model, or adopt prior parameter values to estimate the level of 
contamination. And the wide range of possible distortions for the same correlation-function 
shape resulted in seriously large errors. 

To illustrate how significant these limitations are, and how biased the use of BAO 
measurements of H(z) can be, consider the considerable disparity between Shafer’s 
conclusion that these data strongly favor ACDM over Rh = ct, and the results of the 
Alcock-Paczynski test using only the two most recent measurements that avoid such 
model-dependent biasing (Melia & Lopez-Corredoira 2015). The much more precise 
determination of the Ly-a and quasar auto- and cross-correlation functions at z = 2.34 has 
resulted in the measurement of BAO peak positions to better than ~ 4% accuracy (Delubac 
et al. 2015). A similar accuracy has been achieved through the application of a technique 
of reconstruction to improve the signal/noise measurement of the BAO peak position in the 
anisotropic distribution of SDSS-III/BOSS DR11 galaxies at ( z) = 0.57 (Anderson et al. 
2014). 

LInlike previous measurements, the actual shape of the BAO peak does not affect 
the calculation of its centroid position for these two high-precision cases, both along the 
line-of-sight and in the direction perpendicular to it, when its FWHM is very narrow. The 
peak’s narrowness mitigates the impact of redshift distortions, which affect the peak’s 
amplitude, but not its location. Notice, for instance, that, although the errors for the 
redshift distortions quoted in Table 2 of Delubac et al. (2015) are very large, the relative 
error bars for H(z) are much smaller. In addition, the application of the Alcock-Paczynski 



14 


test to these two precision measurements ensures that any model dependence in the BAO 
data is completely removed because in this test, both the unknown acoustic scale and the 
Hubble constant H 0 completely cancel out, not to mention that this test is also completely 
independent of any possible redshift evolution in the acoustic length. As shown by Mclia 
& Lopez-Corredoira (2015), this test indicates that the concordance model is excluded at a 
99.34% C.L., while the probability that R h = ct is consistent with these data is ~ 96%. 

In other words, the true model-independent BAO data overwhelmingly favor Ry, = ct 
over wCDM/ACDM, in sharp contrast to the conclusions drawn by Shafer (2015). There is 
therefore no merit to his claim that the cosmic-chronometer measurements of H(z), which 
strongly favor = ct over ACDM, should be supplanted by the BAO measurements, which 
he erroneously believes yield the opposite result. 

Shafer (2015) also considers Type la SN data, and though SNe are not the subject 
of the present paper, it is nonetheless useful to point out the bias in this analysis as well. 
These issues have already been discussed and published in, e.g., Wei et al. (2015b), so we 
will provide only a brief summary of the key issues here. 

The statistical analysis of Type la SNe can be improved with the merger of disparate 
sub-samples. For example, the Union2.1 catalog (Kowalski et al. 2008; Suzuki et al. 2012), 
which currently includes ~ 580 SN detections, offers several statistical advantages, but each 
sub-sample comes with its own set of systematic and intrinsic uncertainties. These are 
subsumed into unknown intrinsic dispersions (Tint’s (one for each sub-sample), which makes 
it difficult to fit cosmological models. The commonly followed approach, and apparently 
the one also followed by Shafer (2015), is to minimize an overall y 2 , while constraining the 
y 2 of of each sub-sample to equal unity. It is hardly surprising, then, that the overall yjj of 
is very nearly one. However, a correct statistical approach would estimate the unknown 
(Tint’s simultaneously with all other parameters (Kim 2011; Wei et al. 2015b). For this, 
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the use of maximum likelihood estimation (MLE) has been shown to yield superior results 
(D’Agostini 2005; Kim 2011, though the presence of multiple cr int ’s complicates the analysis 
when the number of merged sub-samples is greater than « 2 — 3. 

The Union2.1 catalog contains at least 17 sub-samples. In reality, therefore, the total 
number of “nuisance” parameters is 20, since all of the <Ji nt ’s need to be recalibrated for 
each independent cosmological model in a truly unbiased test. This is unrealistic. Note that 
this also means the expressions used by Shafer to estimate the AIC and BIC are incorrect. 
Since the (Tint’s themselves are not known a priori, the AIC and BIC must be calculated 
in terms of the likelihood function, not y 2 , as described in Wei et al. (2015b). In such 
cases, it is therefore preferable to work with a few large sub-samples, rather than many 
smaller ones. Fortunately, about half of the Type la SNe in Union2.1 came from the single, 
homogeneous sample known as the SNLS (Guy et al. 2010), and since the same instruments 
and reduction techniques were employed for all 252 of these high-z (0.15 < z < 1.1) events, 
a single cr int characterizes the unknown intrinsic scatter in this homogeneous sample. Notice, 
in particular, that this sample spans the important range of redshifts within which the 
transition from deceleration to acceleration is thought to have occurred. 

This was the approach followed by Wei et al. (2015b), who concluded from the careful 
analysis of the SNLS catalog that the BIC favors = ct over wCDM/ACDM with a 
likelihood of ~ 90% versus only ~ 10%. In contrast, Shafer based his analysis on merged 
sub-samples, with the additional complication of unknown (x int ’s and the use of questionable 
statistical techniques. But he does not adequately explain why his result is completely 
reversed when one uses a single sample with uniform systematics. Surely the outcome of a 
detailed Type la SN analysis, if carried out properly, should be robust enough to emerge 
intact regardless of whether one uses a single large sample, or many smaller sub-samples. 

Indeed, to illustrate this point, we (Melia et al. 2015b) carried out a similar analysis 
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to that of Wei et al. (2015b), this time using both the SNLS and Sloan Digital Sky Survey 
(SDSS-II) events in the Joint Lightcnrve Analysis (JLA) of Betoule et ah (2014). But a 
careful screening of these two sub-samples shows that the SNLS and SDSS-11 SNe have a 
small relative offset in their measured magnitudes, possibly due to slight differences in the 
background subtraction between the two surveys. A systematic shift such as this is not 
adequately handled by introducing sub-sample-specific a int ’s. Instead, a more meaningful 
approach is to use a single sample-wide a int , together with a systematic magnitude offset 
parameter AM 0 ff set between the SNLS and SDSS-II events. The results of this study show 
that the JLA merged SN sample favors = ct over tcCDM/ACDM with a likelihood of 
~ 88% versus only ~ 12%, in complete agreement with the previous study using solely the 
SNLS catalog on its own. 

This is the type of consistency one expects from the analysis of Type la SNe if the 
analysis is being handled correctly and independently of any inherent bias. In particular, 
one cannot use parameters optimized for one model to test another, and the use of a single 
large sample should not completely reverse the results of analysis of more complicated, 
merged catalogs. In summary, Shafer’s (2015) claim, that wCDM/ACDM is preferred over 
i?h = ct by the BAO-measured values of H(z) and the analysis of Type la SNe, is based on 
a combination of improper data selection and a statistically flawed analysis. Instead, the 
results of the present paper show that model-independent measurements of H(z) strongly 
favor R.}, = ct over wCDM/ACDM, in complete agreement with most of the other published 
one-on-one comparative tests between these two cosmologies. 


4. Conclusions 

In reporting the high-z measurements of H(z) using cosmic chronometers, Moresco 
(2015) applied this diagnostic to ACDM in order to demonstrate its power for improving 
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the precision with which the model parameters may be optimized. He found a detectable 
improvement (~ 5%) in the inferred value of ff m and w^ e over previous studies restricted to 
measurements of H(z) at z < 1.75. 

In our previous study using a sample similarly restricted to z < 1.75 (Melia and Maier 
2013), we derived respective Bayesian posterior probabilities of 91.2% for R h = ct and 
8.8% for wCDM/ACDM. In the work reported here, the addition of 4 new measurements 
at z < 0.3 by Zhang et al. (2014) and, especially, the two new measurements at z ~ 2 
by Moresco (2015), have not only confirmed our earlier results, but have strengthened 
the statistical significance of this important one-on-one comparison between competing 
cosmological models. We have found that the model-independent cosmic chronometer data 
prefer the 7% = ct Universe over tcCDM/ACDM with a BIG likelihood of ~ 95% versus 
only ~ 5%. 

We therefore confirm Moresco’s (2015) conclusion that upcoming, additional 
measurements of H(z) with expanded surveys at high redshifts (e.g., Laureijs et al. 2011) 
will provide us with one of the most powerful probes of the cosmic spacetime in the local 
Universe. 

In so doing, we have achieved the primary goal of this work—to independently confirm 
the outcome of previous comparative tests between R h = ct and wCDM/ACDM using 
SNela and BAO data within the critical redshift range (0 < z < 2) where the transition 
from decelerated to accelerated expansion is thought to have occurred. 

These three diagnostics, the SNela Hubble diagram (Melia et al. 2015b), the 
Alcock-Paczynski BAO test (Melia and Lopez-Corredoira 2015), and now the H{z )-redshift 
relation, each probes the cosmic expansion in its own unique way. The fact that all 
three have self-consistently demonstrated that R^ = ct is preferred over wCDM/ACDM 
with comparable statistical significance argues strongly in favor of the zero active mass 



equation-of-state. Our results affirm the influence of dark energy in the cosmic fluid; but 
not in the guise of a cosmological constant. 
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